Abstract 

Kinetics of the first order pfiase transition has been investigated. The 
case when the droplets sizes have a limit value is considered. An analytical 
theory for such process has been constructed and all main characteristics 
of the process have been determined by rather simple approximate an- 
alytical formulas. The limit for a droplets size violates the avalanche 
consumption of the mother phase which is the basic feature of the first 
order phase transition and have been actively used in previous publica- 
tions. Now one can not state that the main quantity of the already con- 
densated substance is contained in the droplets formed under conditions 
unperturbed by nucleation process. This requires to introduce some new 
methods of solutions for kinetic equation which are given in the present 
paper. 
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1 Introduction 

The most interesting feature of the first order phase transition is the global 
kinetics of transformation of a mother phase into a state of objects of a new 
phase. Traditionally the theoretical description of the first order phase transition 
kinetics is given on example of nucleation 0. Then the mother phase is a 
supersaturated vapor and the new phase is a liquid phase. The driving force of 
the first order phase transition is initiated by metastability which is traditionally 
described by the value of a supersaturation (. The last is the ratio of the density 
of surplus vapor and the density of the saturated vapor minus one unit. 

Kinetics of nucleation has been already investigated in the case of the or- 
dinary (unlimited) phase transition. In that case the property of the blow-up 
growth has been actively used . This property means that the rate of growth 
of the molecules number inside the droplet for a given droplet rapidly increases 
in time. It grows so rapidly, that even acceleration of the rate of growth also 
grows in time. This property leads to two important consequences, which allow 
to give an analytical description of the nucleation period. These consequences 
are the following: 

• The nucleation period (i.e. the period of rather intensive appearance of 
supercritical droplets) is well localized in time. 

• The droplets appeared in the first moments of nucleation period contain 
the main quantity of molecules (in comparison with other droplets) during 
the nucleation period. These first moments can be called as the pre- 
nucleation period. During the pre-nucleation period the droplets don't 
essentially perturb the state of the system. This allows to calculate the 
characteristics of droplets appeared during the pre-nucleation period. 

For example, the regime of the unlimited free molecular growth has the men- 
tioned feature. Namely here the number of molecules v inside the supercritical 
embryo (i.e. inside the droplet) grows in time t as 

diiyy/^/dt^C/T, a = 3, (1) 
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where r is a characteristic constant which has a sense of a mean time between 
colhsions of a given molecule in saturated vapor with other molecules. This type 
of growth is wide spread in nature. Later we shall use this type of growth as a 
model for concrete formulas. 

In previous theories for nucleation kinetics (see rewiev in pQ) it was sup- 
posed that the sizes of droplets can become in principle infinitely big. Here we 
consider the opposite situation. In many systems the process of the droplets 
growth is limited by some artificial restrictions. The number v of molecules 
inside the droplet can not be greater than some value vum, which is an external 
parameter of the theory. We imply that strongly exceeds the number in 
the critical embryo under the supersaturation which corresponds to relatively 
intensive appearance of droplets (this condition will be specified) . The existence 
of a limit size takes place in the processes of cementation, of condensation in an 
external lattice, morphological transitions, etc. This type of transition is rather 
wide spread in nature. Kinetics of this process will be investigated in this paper. 

The limit of droplets sizes leads to the loss of the blow up property of growth 
in the global kinetics of nucleation. From the first point of view it seems that 
the introduction of a limit size can not essentially change the approaches of 
description, but it is an illusion and below some principally new methods of 
solution of kinetic equations have to be presented. 

This publication continues the set of articles devoted to the construction 
of the global kinetics of a phase transition (see PP). To grasp the current ap- 
proach it would be useful to come through these papers. We use all physical 
assumptions from these papers except explicitly mentioned ones. All definitions 
of physical characteristics can be taken also from ^ . 

Traditionally two characteristic types of external conditions are considered. 
The first type is a decay of a metastable state [21 . Here at the initial moment 
of time we have the metastable phase with no objects of a new phase in the 
system. Later there will be no external influence on the system. The process of 
nucleation begins and the appearance and growth of droplets cause the fall of 
supersaturation which stops the further intensive nucleation. 

The second type is nucleation under the smooth variation of external con- 
ditions. At first moments of time the increase of supersaturation by the action 
of external conditions is greater than the decrease due to the mother phase 
consumption by the droplets. But the vapor consumption has a very rapidly 
increasing intensity and later the mentioned decrease will overtake the men- 
tioned increase and the supersaturation will go down. The process of intensive 
nucleation will stop. We are going to describe kinetics of the process in both 
situations. 

We shall use the following law of droplets growth for p = v^l'^: 

'^P ri t / -1/3 
= C/T for p < pH„, = Vil^ 

dp n f ^ - 1/3 
— = for p> pum = ^iL 
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Rigorously speaking it would necessary to determine the variable which has 
the velocity of growth independent on the value of this variable. We shall act in 
another manner: at first we suppose that there is no limitation on growth and 
then calculate the number of molecules in a droplet by a "curved" formula. 

Here it is possible to introduce the following variables: The variable y is the 
number of molecules in the droplet which is growing without any limit by an 
initial law . The variable z is the value of y for droplet which was born at 
"initial" moment of time. Certainly, for every droplet v < z. The variable x is 
defined as y = z — x. Then we see the correspondence between t and z and the 
correspondence between t and x. The values z and x are now equivalent. 

The number of molecules v inside the droplet has to be calculated now as 

V = {z - xf for z-x < pum 

and 

V = Vlim for Z - X> piim 

We see that here there is no unique analytic dependence and this will seriously 
complicate the evolution equation. 

We have also to mention that the property of the effective size of growth 
which was observed in recent publication in construction of kinetics for the 
law ^ with small a < 1 also doesn't take place. For droplets with^ v ^ 
such a property is observed, but for droplets with v < vn^n is the opposite 
property of the blow-up regime takes place. Namely the situation where both 
regions are important in substance balance is interesting. 

All analytical constructions are based on the following approximation for the 
rate of nucleation / as a function of supersaturation: 

/(C) = /(Co)exp(-p(C-Co)) (2) 
Here Co is the base of decomposition and the positive parameter p is defined as 

dF, 

where is the free energy of a critical embryo. The last approximation is not 
more than the Klapeiron-Klausius approximation, but it is written for special 
variables. The validity of this approximation has been discussed in 

2 Decay of metastable state 

Now we shall consider the situation of decay. The initial value of supersaturation 
will be marked by The base of decomposition Co for function C, can be chosen 
as $. 

^Here y is imaginary value calculated purely according to tTl . 
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We skip the derivation of kinetic equation because it is quite similar to the 
ordinary case without the limit size and can be seen in ^ or in |Sj . The kinetic 
equation for the evolution can be written in the following form: 



$ = C(^;) + 0(2 - zum){A / (z - xy exp(p(C(a;) - <i>))dx + 

(3) 

Avum / exp(p(C(a;) - + Q{zum - z)A {z - x)^ exp{p{({x) - <i>))dx 

Jo Jo 

Here A is the amplitude value of spectrum, is the Heavisaid's function, zum = 
piim- Solution of this equation will give the behavior C as a function of x. It 
allows to calculate the number of droplets (in renormalized units) 



N{z) = / exp{p{Cix) - <^>))dx 
Jo 

appeared until the moment t{z) (i.e. until the moment when the imaginary front 
side of the spectrum attains the value z). The most interesting value is the total 
number of droplets N{oo) and the error in this value will be the measure of the 
accuracy of approximate description. 

In contrast to the decay with the unlimited growth we shall have here the 
several stages of the process of decay and the long continuous tail. The meaning 
of these parts of the process will be clear later and now we shall start with 
investigation of the first period and have to put initial conditions. Namely, 
t — 0, C — ^ and there is no droplets in the system. 

Now we shall introduce the value of natural length Aix as the characteristic 
width of spectrum without the limit of growth. For this situation we have the 
following kinetic equation 

$ = C{z) + A f {z-xf exp(p(C(a;) - ^))dx 
Jo 

which results in the following characteristic width of the size spectrum 

/ - exp(p(C(a;) - $)) 

namely 

One can see that Aiz is a function of $ since p = p{^), A = A{^). It is 
possible to see three situations: 

(A) AiZ < ZHm, 

(B) Aiz > ZH. 

(C) Aiz « zii- 



•Lira ; 
lim • 
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2.1 The case Ai2 <ti 

In this case kinetic equation |j2Jl can be reduced to 



G{z) {z-xy eyip{~pG{x))dx 

Jo 

for the function ^ 

G^A I {z-x)^ exp(p(C - <i>))dx 
Jo 

which has the sense of the number of molecules inside the new phase. 

This equation is known from unlimited case and can be solved by iterations: 

Go -0 

G{i+i) (z) = A (z - x)^ exp(-pG'(i) {x))dx 
Jo 

Calculation gives G(i) = Az'^/A and already the second iteration gives the prac- 
tically accurate number of droplets 



i-x'')dx = 0.9 



7V(2)M = y^ eM-pGi{x))dx^{—y/''B, B = J^ exp( 
The relative error 

|iVM~jV(2)MI 

^ ^ iV(oo) 

is less than 0.05 in every case which can be analytically proven. 
The high accuracy is also attained also for 

3! f°° 
Li ^ 7777. TT^ / exp{-pG)dx 

calculated on the base of the first iteration i.e for 

3! f°° 
Lt « L, (2) = — -tA / a;*exp(-pGi)da; 

These values will be important later. 

But this doesn't complete the investigation as in the case of unlimited 
growth. Earlier or later the droplets sizes attain h'um ^^nd the droplets stop 
to grow. Then the final value of the supersaturation will be 

C/m 1 = N{oo)viim 

Here index 1 points that we investigate the first part of the nucleation period. 

The time of attaining can be easily calculated on the base of monodis- 
perse approximation^. The last approximation leads to equation 

dz 

AN{oo)z^ = $ - (5) 



The validity of monodisperce approximation can be proved analytically. 
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which can be easily integrated (as the first order differential equation without 
explicit dependence on the argument): 

r dx _ ^ 

Jo ^ ~ AN{oo)x^ ~ 

The integral can be taken analytically and gives z{t). Condition 

gives the value of the end tjin of the first part of nucleation. One can compare 
tfin with the duration of intensive nucleation during the first part, i.e. with 

Alt = Ai2t/$ 

When the inequality Aiz ^ z/i„ isn't too strong then the equation Q isn't 
suitable because the monodisperse approximation lying in the base of (O isn't 
valid even at z w zjin- Then one has to use instead of (jSJ the following equation 

^L,z^ = <i>-T— , z|t=o = (6) 

i=0 

The methods of solution of this equation are quite similar to the methods of 
solution of © (they are presented in P). 

While approaching to zji„ one has to take into account that Li begin to 
change because integration can not be done out of z < Zfin- These changes can 
be taken into account by the perturbation technique (because the main influence 
is given by Lq which remains unperturbed). 

Several important statements can be proved here. They are: 

• (A) The process of nucleation takes place in the quasistationary conditions. 
The time ti of essential perturbation of the droplets distribution, which 
can be estimated as Azit/$ strongly exceeds the time tg of establishing 
the stationary state in the near critical region (more accurate in the region^ 
y < (2^3)i/c). 

• (B) The main role in the vapor consumption is played by the strongly 
supercritical droplets, i.e. by the droplets with i' > {2 ^ 3)i'c. These 
droplets grow regularly with the mentioned law of growth. 

• (C) The process of intensive nucleation takes place at $ > C > $ — 1/p. 

• (D) Later the new periods of nucleation (may be not so intensive but long) 
will take place. 

The statements (A), (B) are necessary for construction of the kinetic equa- 
tion. The statement (C) allows to use the approximation The statement 

^To calculate ts see |7|. 
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(D) requires to continue the consideration of the nucleation process which will 
be done below. 

The last step to do is to give approximation for the number of droplets. 
After the evident rescaling the kinetic equation can be presented as 

G = 4 [ {z-xf exp{-G{x))dx 
Jo 

The number of droplets appeared until z can be calculated as 

N{z)= [ exp{-x^)dx 
Jo 

We shall use the following approximation 

Nappriz) = 9(0.91 - z)z + e(z - 0.91)0.91 

The relative error \Nappr — N\/N as function of x is drawn in Figure 1. We see 
that it is small. 

2.2 The case Aiz ^ 

This case is much more simple than the previous one. Here the kinetic equation 
can be written as 

$ = C{z) + Aviim / exp(p(C(a;) - ^))dx 

Jo 

or 

G = Aviim / exp{-pG{x))dx 
Jo 

for 

G = Auiim / exp(p(C(x) - 

Jo 

It can be solved very easy. Namely, after the differentiation we come to 



dG 

— = Aviim exp(-pG(2;)) 
which can be easily integrated 

i(exp(pG(z)) - 1) = AuiimZ 
P 

The spectrum has the form 

f ^ Aexp{-pG{z)) 
One can see that here the number of total droplets 



N{oo) = / f{x)dx 
Jo 
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is infinite. It is certainly a wrong result. The reason is inapplicability of ap- 
proximation (0) • The head of the spectrum is described quite satisfactory, but 
the tail isn't described well. Really at n = rioo the nucleation according this 
approximation doesn't stop. To describe the tail we have to use the precise 
equation for the stationary nucleation rate 

= Zexp(-F,) 

where Z is the Zeldovitch factor [HI, Fc is the free energy of the critical em- 
bryo, taken in some approach but without approximate decomposition Q . For 
example, in capillary approximation 



2a 



3 = ' ^ 31n(C + l) 

where a is renormalized surface tension. Then 

ft 

dm 



or 

The last equation can be easily integrated 



dt 



= t 



l^limlsiCit)) 
^ 6^P( 271n"fc+l) )^^ _^ 

One can get an analytical solution for the tail of spectrum having noticed 
that at small C one can decompose ln(C -I- 1) ^. The values of ^ can not be 
too small, at least fdC) < t^Um and the behavior of Z isn't singular. So we can 
reduce the last equation to 



B2 J cxp(B3C"')rfC = t + Bi 



with known constants Bi, i = 1,2, 3. The integral can be reduced to the error 
function. This solves the problem of explicit solution in terms of standard special 
functions. 

Here one has no need to observe the property (B) because this property 

1 /3 

begin to be violated only when i/^^^ strongly exceeds Azi and then practically 
all droplets attain the final values. As for the quasistationarity (property (A)) 
we see that \d(/dt\ decreases in time and if (A) is observed in the initial moment 
of time it will be observed also later. The absence of the property (C) was taken 
into account just above. Instead of (D) one can state that this stage will be the 
last stage of the whole nucleation process in the system. 
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2.3 The case Aiz zum 

The intermediate case is the most complex one and has to be investigated on 
the base of the already presented solutions. Let the kinetic equation be rescaled 
to have 

G = 

Q{z - -2iim){4 / {z- x)^ exp{-G{x))dx + Ai^nm / exp{-G{x))dx}+ 

Jz — Zlim Jo 

Q{ziim - z)A {z- x)^ exp{-G{x))dx 



Certainly, here zum will be also rescaled. 

Now we shall give the approximate method to solve this equation. The value 
G can be presented as 

G = g++g- 

where g+ is the rescaled number of molecules in the droplets which have already 

attained , g- is the rescaled number of molecules in all other droplets. 

For further constructions we have to write approximations for (/_ and 
We shall start with g- 

The function g_ can be easily found on the base of iteration procedure 
presented for the case diz <^ zum- Then 

POO pZ 

g- = Q{z-/^-i_z)AA I {z-xf ex:p{-g--g+)dx+Q{Aiz-z)AA / {z-xf ex.^{-g--g+)ddL 
Jo Jo 

Here A is the amplitude which can differ from 1 due to the possible big value 
of g+. At the current moment this remark isn't too clear, but in any case we 
can say that for further purposes it is convenient to conserve A here. 

Now we shall explain the last relation. 

At ^; < ziim we have g+ = and 



g-{z) = AA / {z — x)^exp{—g-{x))dx 
Jo 

It is important that the last equation is the closed one and it is identical to 
the case of unlimited growth. So, it can be successfully solved by iterations. 
Particularly, (?-(o) = 0,5_(i) = Az'^, etc. 
At ^; > ziim we have 



AA / {z — x)'^ eyip{—g-{x) — g+{x))dx 

J Z — Zli^. 



Now we give the qualitative picture of phenomena. At the first moments of 

time the spectrum in the region p < pum is being formed and g+ = but g_ is 
growing rather rapidly ~ z'^. Later the substance is going to be accumulated in 
droplets with p = pum- The growth of g+ occurs. 
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If the size of cut-ofF z^m is essentially smaller than Aiz = A~^^^ then g_ 
accumulates not so many molecules of substance g- = g-{in) = ^-^fim ^ ^ ^^'i 
later due to the decrease of nucleation rate the value of g- will also decrease, 
i.e. g- < g-{in) In this case 3+ grows as 

(•2-ZHm 

G+ w Z(im4A / exp(-g+(a;) - g-{x))dx 

Jo 

It will grow slower than 

9+{in) = ZUm^A{z - ZUm) 

This allows to state that 

dz 

Then the variation Sg+ in 5+ during "the time" of establishing of the quasista- 
tionary state in the region [0, zum] can be estimated as 

5g+ < ^ziim = ^Azt„^ < 1 
dz 

In the expression for g- for z > zum the last inequality allows to take exp(— 
away from the integral. 
Then approximately 

g-{z) = 4Aexp(-g+(z)) / {z - xf exp{-g_)dx 

•I Z — ZHm 

and on the base of initial iteration approximation 

g-{z) = Aexp(-fi(+(^))4m 

Now we shall get expression for g+. We shall get for g+ the following ap- 
proximate equation 

g+ = 4:Azfi^ / exp(-5_(a;) - g+{x))dx 

Jo 

Here z is the current " moment of time" , zum is the " moment" when the spectrum 
attains zum- 

Here we are interested in the other characteristic scales of time which are 

greater than l/{4,Aziim). At these times we can suppose that exp(— g_) is 
approximately constant and take it out of the integral. Then approximately 

fif+ = 4A4™exp(-5_) / Qxp{-g+{x))dx 
Jo 

When g_ is essential (i.e. exp(— (7_) is small) then isn't necessary and we can 
take out exp(— But when g- is small, then exp{—g-) = 1 and we can also 
take it out. The substance consumption in g- occurs in the avalanche manner. 
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Wc can get solution of this equation quite analogously to the case Aiz ^ zum 
but here the initial conditions will be another 

g-\x=Z-ZUm = 

Solution will be 

^ = 4A4m exp(-5(-) exp(-5+) 

and 

exp(5+) - 1 = 4^4m exp(-3_)(z - zum) 
The spectrum will be 

/ = 4^exp(-g_)/exp(g+) 

or 

/ = 4Aexp(-5_)/(l + 4^4„ exp(-fif_)(2; - znm)) 
In the initial approximation we consider A = l and get 

/ = /o = 4exp(-4m)/(l + exp{-g_){z - zum)) 
Then we take g- = zf-^ and get 

/ = /o = 4 exp(-4m)/ (1 + 44m exp(-4m)(^ - Zlim)) 

It is necessary refine expression for the spectrum obtained in the zero ap- 
proximation. We have to take into account that is formed under the influence 
of g-. But the spectrum in the last expression is written through £f_ (but not 
through Here is the difficulty. We shall come it over via the special ap- 

proximate approach. At first we assume that g- (but not the g-^-) is formed 
with the renormalized intensity. We take exp(— a^m) instead of ex.p{—zii^i). 
It is evident that we do this only in the first factor which is responsible for the 
influence of g^. So, wc get 

/ = 4exp(-a4m)/(l + 44m exp(-4m)(2 - ^Hm)) 

For a we have to take the renormalized amplitude. The natural candidate for 
this amplitude is /q. As the result we come to the following expression for 
spectrum 

/i = F(/o) = foK 

where 

exp(-4m) 

^'^1H--/1,J 

° 1 + 44m exp(-4m)(^ - Zlim) 

Now let us realize what have we done. In reality wc had to take into account 
the inverse influence, i.e. to take into account that <?+ is formed under the 
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influence of g-, but not the influence of g+ on To take into account the 
inverse influence we shall use the inverse transformation . So, approximately 

We can calculate this transformation rather easy 

This will solve the problem of the adequate construction of the approximate 
spectrum. 

The error of this approximation in the total droplets number is drawn in 
Figure 2. Certainly, we can not calculate here the total number of droplets (it 
will be inflnite and has to be corrected by the way described in the previous 
subsection). As the final value Zf in the magnitude 10 has been chosen. The 
last belongs to the asymptotic period. One can see that the relative error in the 
droplets number is small. 

Statements (A), (B), (C) take place here. Statement (D) isn't necessary. 

2.4 Further nucleation 

The case Aiz ^ requires the analysis of the further nucleation. This 
analysis isn't too complex. On the base of C,fin i we can calculate the length 
A22: as 

A2Z ^ Aiz(C/m 1) 

i.e. we use the same formula but with C,i f m instead of In expression there 
is an amplitude of spectrum A which is a very sharp function of supersaturation. 
So, we see that 

Aiz < A22; 

This strong inequality allows to neglect the time interval [0, Aizr/<i>] in the 
nucleation at C2- Then the situation is reduced to the nucleation with initial 
supersaturation (,2- We have to repeat the steps of the three previous sections. 
Again we can come to the possibilities (B) and (C). The possibility (A) cannot 
take place because the inequality 

exp(iVia;fj„/4) < xum 

doesn't take place. Here the l.h.s. is A2Z in rescaled units, A^i has to be put 
to 0.9, the coefficient 4 in denominator came from the power 1/4 in Q), the 
parameter xum has to satisfy xum > 2 3. 

So, we see that the structure of spectrum isn't too complex - the plateau, 
the region of partial collapse of surplus substance and may be the tail. 

The analysis of nucleation after instantaneous creation of a metastable state 
has been completed. 
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3 Nucleation under the smooth behavior of ex- 
ternal conditions 



The straightforward generahzation of kinetic equation to the case of smooth 
change of external conditions leads to the following kinetic equation 

$(z) = C(^) + (7) 
A {z- xf exp(p(C(x) - <^^))dx + Avum / exp(p(C(a;) - 



Here ^[z) is ideal supersaturation, i.e. the supersaturation which should take 
place in the system without any processes of vapor consumption and heat release 
effects. The value of $ at some characteristic moment will be marked as 
We shall choose later. The number of droplets (in renormalized units) will 
be calculated as 



N{z) = / exp(p(C(a;) - ^))dx 

J —oo 

Under the unlimited growth of droplets the period of nucleation is well lo- 
calized in time. Moreover it is rather short in time which allows to linearize the 
ideal supersaturation during this period. Here the process of nucleation cannot 
be localized in time, it is seen simply by the fact that the number of droplets in a 
liquid phase is limited from above by Nvum ■ Then to compensate the " action" 
of external conditions it is necessary to have nucleation again and again. 

To present concrete calculations we suppose that the ideal supersaturation 
can be linearized 

$ z = + -^U=oz 
dz 

but now is is no more than a model. Later we shall see how to construct the 
theory for rather arbitrary behavior $(2:). 

But even with this linearization the situation here can not be solved by a 
simple generalization of the ordinary iteration procedure because here process 
isn't limited in time. 

3.1 Rescaling 

Kinetic equation (UJ can be written in the following form 

/z-xii„-i l-Z 
XH.mf{x)dx+ / {z-xff{x)dx 

with xum = Zlim and a spectrum f{x) 

f{x) — aex.p{bx — g{x)) 

It is essential that after the evident rescaling one can put a = 1 and 6=1. 
The third parameter xum remains and this leads to the absence of universality. 
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Restriction a = 1, 6 = 1 doesn't lead to the special meaning of the point 
z = (in homogeneous condensation z = was the point of the maximum of 
supersaturation, in heterogeneous condensation z = was the point where the 
half of droplets have been appeared). One has simply require that <i>(z = 0) isn't 
too far from the maximum of supersaturation. This fact can be also analytically 
proved here. 

When xiirn = oo as one can get the half-width of spectrum by condition 

Az = 1 

3.2 Pulse regime 

When xum 3> 1 the process of nucleation can be described rather simple. To 
describe the first peak of the nucleation intensity we have to solve equation 
without limit of growth. This case is well known and can be solved both by 
iterations [S] and by universal solution (j. Thus, we get Ni - the number of 
droplets formed in the first peak. One can see A'^i ~ 1. Then the supersaturation 
^ continues to fall until z xum ■ The inequality xum ^ 1 allows here to speak 
about the well formed peak of nucleation intensity and to use the standard 
procedure of unlimited growth. At z ^ xum the supesaturation begins to grow. 
At z ~ xiim we have / ~ exp(a;Hm — ^i^um) ^1- At z ~ Nixf^^ the second 
peak of nucleation begins. 

The second peak of nucleation can be described absolutely analogously to 
the first one. So, we have the set of identical peaks. The nucleation occurs in 
the pulse regime. 

This regime is shown in Figure 3. 

3.3 Smooth regime 

The opposite situation when znm <C 1 can be also solved analytically. Here 
kinetic equation can be written as 

exp{p{({x) - (8) 

-oo 

and ^ 

$(z) = C{z) + Avu^ / exp(p(C(x) - (9) 

J —oo 

After differentiation with account of linearization we have 

^ = Avuyn exp(6z - _pg(z)) 
dz 

In our units p = 1, Avum = 1,6 = 1. Integration with initial conditions 

g[z = -co) = 

gives 

exp(g(z)) = Aviim exp(z) + 1 
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This value will be marked as Qst- Then the spectrum / will be 

f = fst = exp(a; - ln(exp(ar) + 1)) (10) 

At a; — »• 00 we see fst um = 1- 

This spectrum f^t and the example of numerical solution are drawn in Figure 

4. 

These limit situations will form the base for description of nucleation in all 
situations. But here we have to stress two features: 

• The limit solutions have to be radically modified. The force of the al- 
ready presented methods for these solutions isn't sufficient to give the 
total description. 

• Description of the base of the iteration method can not be suitable because 
the process isn't restricted in time. 

At first we shall present the general structure of methods and then we shall 
show how to determine the elements included in these approaches. 

3.4 Advanced pulse regime 

To describe situations with xum > 2 and \xiim — 2| <C 2 one has to use the 
rescaled peaks approach. This approach is based on the following simple con- 
siderations. 

Due to the previous analysis we shall approximate the spectrum by several 
peaks. Every peak is formed by some external effective source (initiated both 
by external conditions and by the action of the tail of the previous peak). The 
distance between peaks of nucleation is approximately 

Afx = Nxl^ 

where N is the number of droplets formed in the previous peak. Certainly, 
between different peaks the distance can be different because N will be different, 
but the suitable approximation is to consider N to be the ideal value, i.e. at the 
first peak. The last formula will be rather evident if wc notice that the action 
of external conditions has simply to compensate the loss of substance which is 
equal to Nxl^. 

Here we shall present a recurrent procedure to calculate peaks of nucleation. 
Suppose that we have described the current peak and know the coordinate of 
maximum Xm- 

Now wc have to rcscalc the amplitude of the next peak. Wc neglect the 
non-linear behavior of effective source. Then the next peak will be similar after 
rescaling to the previous one. But without rescaling in the same units the next 
peak will have a new (ordinary smaller) height and another width (ordinary 
wider). So, now we have to determine the new intensity h' of the effective 
source (the initial intensity was 6o = !)• 
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We know the expression for intensity of vapor consumption by droplets i.e. 

dg/dz: 

^ = 3A f {z — exp(x — g{x))dx 

The subintegral function has the form 

{z — x)^ exp(a; — g{x)) 

and approximately has maximum a.t {z — x) = pa = 2 for xum > 2 and {z — x) = 
Pa = xiirn for xum < 2 (hcrc wc can put 2 = at maximum because this choice 
will result only in parameters but not in the spectrum form). 

The current peak have coordinate Xm, the next peak will have coordinate 
Xm + ^fX. The maximum of subintegral expression for d,p/dz for the next peak 
will be attained at Xb = Xm + ^fX — aJo.The main effect is connected with a 
property that a.t x = Xb the essential quantity of droplets from the previous 
peak has attained the value xnm- B\it not all droplets have attained this value. 
So, the supersaturation ( lies higher than the supposed value 



602: - Nxl 



3 

im 



which would ensure the similarity of the next peak to the previous one. 

The supersaturation at the corresponding coordinate Xc = Xm — Xa for the 
previous peak was lower than at the current peak. 

We suppose that at x = Xb + Xa = Xm + Aj-x approximately corresponding to 
the coordinate of the next peak all droplets of the previous peak will attain the 
limit value Nxf^^ and there will be no difference between the value of effective 
ideal supersaturation (the supersaturation imaginary formed without droplets 
of the last peak taken into account) fl at Xm and at Xm + AfX 

^{xm) = ^{xm + A/a;) 
The effective intensity b' of the external source can be approximated as 

i}{Xm + AfX) - Ct{Xm + AfX - Xa) 

~ 

Xa 

and have to be compared with 

^{Xm) ^{Xm Xa) 



^0 

Xa 

Then taking into account 

\Q.{Xm + AfX - Xm) - ^{Xm + Afx)\ <C 1 

and 

|^(*^m) ^{Xm ^{Xm Xa)\ Xa ^ 1 
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one can come to 

Xa 

where 

AC = Q{Xrn + A/a; - Xa) - ^{Xrn ~ Xa) 

Since bo ~ I one can come to 

b'^l-^ 

We can express AC as 

AC = ln/(xb)-ln/(x,)=ln4[^ 

f{Xc) 

where / is the spectrum corresponding to the external ideal supersaturation. 
Since Xa > 1 one can see that the effect of the influence of the current peak is 
negligible and 

f{xb) = f{xb) 

f{Xc) = f{Xc) 

So, finally we have the following formula for 6': 

b' = 6o(l - -)b 

Xa 

where Ca is the difference of supersaturations at a; = a;c = Xm — Xa and at 
X = Xb = Xm — Xa + ^fX Calculated for the previous peak. Here Xm is the 
coordinate of the maximum of the supersaturation. The difference Xb — Xc can 
be calculated as 

Xb-Xc^ ln( {^^''| ) 

fiXc) 

This completes the procedure. The first peak calculations have to be per- 
formed in explicit manner. 

The calculations are presented in Figure 5. We see that the calculated peaks 
are very close to the real solution. 

Even such simple rescaling brings a rather accurate result. When xum > 2 
the accuracy will be better. It can be proved analytically. The simple account 
of variations of AfX will give the better accuracy also. All other errors can be 
taken into account by the standard perturbation technique. 



3.5 Advanced smooth regime 

We see that already at xum — 2 the first peak lies not so far from the stationary 
solution H10() . Later all peaks will be even smaller than the first one. So, we can 
consider (/ — fst)/ fst as the small parameter and linearize the kinetic equation 
over this parameter. Then this solution can be solved analytically and will give 
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the oscihations relaxing to fst- We can avoid these long analytical derivations 
and simply draw the relaxation oscillations 



OSC — 



fst + ki exp{{x - Xra)/k2) COs((x - Xra)/^fX) 



with two parameters ki and k2- The arguments for derivation of Ay a; takes 
place here also. 

Parameters ki and k2 can be determined by a coincidence of approximate and 
real solution at maxima of two first peaks. We choose maxima for coincidence 
of real solution and approximate solution because namely here the intensity of 
nucleation has maxima. The first local minima demonstrates the deviation of 
approximate solution from the real one but here the intensity of nucleation is 
small. 

This solution is drawn in Figure 6. One can see that precise solution is very 
close to approximate one. 

3.6 Iterations to determine parameters 

Now we know the methods to describe nucleation adequately, but parameters 
in these approximations are unknown and our next task will be to determine 
these parameters. To do this we have to know the evolution during several 
first peaks of nucleation. In the advanced pulse method we have to know the 
evolution during the first peak including the back side of the peak and in the 
advanced relaxation method we have to know the positions and values of two 
first maximums. These characteristics will be determined on the base of iteration 
methods. But here the iteration methods can not be got by direct generalization 
of already presented procedures in the case of condensation under the pure free 
molecular regime of vapor consumption [Hj, 0). Really, even in the pure free 
molecular regime it is difficult to describe even the back side of spectrum on the 
base of standard iterations and here we need to describe the second peak. The 
limitation of growth will also diminish the converging power of iterations. So 
we need to reexamine the iteration procedure. 
The evolution equation can be rewritten as 



.9 = 



•Z 




@{Z - Zlim){ 



{z — x)^ exp(a: — g{x))dx + vi 



exp(a: — g{x))dx}+ 



Qiziim - z) / {z - x)^ exp{x - g{x))dx 
Jo 

The initial standard iteration procedure can be determined as 





gi+i = Q{z-ziijn){ {z~x)^ exp{x~gi{x))dx+i^H„i / 



ex.p{x—gi{x))dx) + 
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It will converge at every initial approximation. 
As suitable initial approximations one can propose 



50 = 



or 



= 9st 



at X < ln(/s4 u„i), gst = x at x> \n{fst Um)- 



Since here it will be impossible to take only two first iterations there is no 
special significance what approximation we shall use. We prefer to use the first 
one because here the chain of inequalities 



for every fixed x is observed. These iterations arc shown in Figure 7. 
Decomposition in iterations 

The analytical calculation of iterations will stop at the second step. For 
the unlimited free molecular regime this number of iteration is sufficient, but 
here we have to calculate the next iterations. It is possible analytically if wc 
choose the decomposition of exponent in the subintegral function. But here the 
following difficulty appear. The value of g goes to infinity and it isn't possible to 
decompose exp(— £f(a;)) over g{x). The other possibility is to decompose exp(a; — 
g{x)) over x — g near some characteristic value of the exponent argument. But 
at a; ^ —00 the value of g goes to zero and the spectrum looks like exp(— 
Then it is impossible to use the second decomposition at initial moments of time 
(i.e. at negative x). But this very period drives the evolution during the first 
peak of nucleation and, thus, plays the main role in nucleation description. So, 
the pure second approach can not lead to a suitable result. So, we need to give 
here more sophisticated approach. 

We shall act in the following manner. Already the first iteration gives us 
the approximate position of the first maximum Xmax i and the value of this 
maximum fmax i- Then we propose the following procedure 

• When X < Xmax 1 then 



Here the first four terms are taken into account. It is sufficient to ensure 
the high relative accuracy. 

• When X > Xmax 1 then 



90 <g2< ••• < 921 < ... < g < ... < 92i+i < ... < 93 <9i 



(11) 




4 



4 



(a: - g - Cbaae)' 



exp(a; - 5) « exp(C6ase) ^ 
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where 



(base = ln( 



max 



1 + fst lim 



) 



2 



is the base of decompositions. 

Here decomposition of the corresponding exponent is Umited also by the 
first fom terms. 

Then the iterations can be calculated analytically at every step and they give 
the adequate approximation for the solution. They converge to the solution of 
kinetic equation where instead of cxp one should write the presented approx- 
imation. Already the four first terms of decomposition of the corresponding 
exponent ensure the high accuracy. 

These iterations are shown in Figure 8. 

We see that already the several first iterations give the practically precise 
positions of the first maximum and the first minimum fmin i at Xmin i- Later 
we see that iterations converge but every new iteration gives the " step" of coin- 
cidence with the real solution which becomes smaller and smaller. The reason 
of this failure is the following: really, the precise solution doesn't go far from 
fst lim and decomposition works, but every iteration inevitably goes to zero or 
to infinity and on every iteration the decomposition does't work. This error 
leads to crisis in convergence of iterations. So, we need to modify iterations. 

Cut-ofF of iterations 

Now we know fmax i and fmin 1 and can use these characteristics. One can 
analytically prove that for every x > Xmax i 



Then we can require that before we use decomposition the following cut-off 

• When exp(a; — gi{x)) > fmax i we should take instead of exp(x — gi{x)) 
the value the constant fmax i and fulfill no further decompositions'*. 

• When exp(a; — gi{x)) < fmin i we should take instead of exp(a; — gi{x)) 
the value the constant i and fulfill no further decompositions. 

• When fmin 1 < exp(a; — gi{x)) < fmax i no special actions are required 
and we have to make the mentioned decompositions. 

The convergence of new iterations is higher than the convergence of previous 
iterations. New iterations evidently converge to the real solution^. Also it is 
very easy to estimate rigorously the error of decomposition because now we have 
an estimate fmin i < fi < fmax i for every x and for every i. 

To complete the description of iterations one can prove the following impor- 
tant statements 

^One can see that iterations escape from the real solution, 
^with corresponding decomposition. 



fmin 1 ^ ./ ^ fi 



max 1 
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• The cut-off iterations without decompositions converge to the real solution 
of kinetic equation and they converge faster than the first-type iterations. 
The chains of inequalities Hll|l remain valid for the cut-off iterations with- 
out decompositions. 

• The iterations with decompositions converge to the solution of kinetic 
equation where function cxp is treated according to the mentioned de- 
compositions. The chains of inequalities remain valid (here g is so- 
lution of kinetic equation where the function exp is treated according to 
the mentioned decompositions)^. 

For practical needs one can show that for xum < 2 already the sixth iteration 
ensures the correct values of fmax i , Xmax 2 and fmax2 1 Xmax 2 (the value and the 
coordinate of the second peak)^. This is all we need to construct the relaxation 
oscillations. 

Behavior of the cut-off iterations is shown in Figure 9. One can see the high 
rate of convergence. 

When xiim > 2 we can use the advanced pulse method and here we need to 
describe only the first peak. It can be done by the first four iterations (Really, 
if it would be necessary to get the coordinate of the second peak we might have 
difficulties because of the very small value of fmin i)- 

This solves the problem to get the parameters of approximation and com- 
pletes the description of the nucleation under the break of the blow-up growth. 

^To show this we have to see that decomposition of exp conserves the monotonious prop- 
erties. It can be shown if we differentiate the approximation and note that the derivative of 
approximation is the approximation for the derivative of exponent. Then with the necessary 
terms taken into account and with account that fmin 1 < /i < fmax 1 this derivative can be 
made positive. Then the monotonious properties remain even after decompositions. 

'^When xiim < 2 then the iterations will approximate these values even better then at 
Xiim = 2. It can be proved analytically. 
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Figure 1 

Relative error e of initial approximation 
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Figure 2 

Relative error e of the special approximation 
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Figure 3 

Peaks of spectrum in impulse regime xum = 2.6 (if it will be greater we 
can not sec the difference between the precise and the approximate solution). 
The variable z is going until 30. Approximate peaks lie higher than the precise 
solution, The last peak lies at greater z than the precise solution. 
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Figure 4 

Peaks of spectrum in smooth regime xum = 1-25 (if it will be smaller we can 

not sec the difference between precise and approximate solution). The value z 
is going until 20. The precise solution oscilates, the approximate solution is flat. 
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Figure 5 

Peaks of spectrum and advanced peaks at XHm = 2.0. 

The value z is going until 25. Approximate peaks are separate lines. The 
global continuous line is the precise solution. 
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Figure 6 

Spectrum and the relaxaiton solution at xum = 2. 

The vahic z is going until 25. Approximate solution has more symmetric 
peaks than the precise one. Only in the first negative gap one can see the 
difference between approximate and precise solutions. 
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Figure 7 

Precise solutions and iterations (calculated numerically) without decompo- 
sitions for xiim = 2 (the worst situation). Iterations arc going near precise 
solution and then fall very rapidly down or grow up to infinity. The higher is 
the iteration number the longer it is going near the real solution 
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Figure 8 

Precise solutions and iterations with decompositions (calculated analyti- 
cally). Here xiim = 2. (the worst situation) Iterations arc going near precise 
solution and then very rapidly fall down or grow up to infinity. With the growth 
of iteration number there is practicUy no effect in the lenght of coincidence be- 
tween the iteration and the precise solution. 



30 



/ 


\ \ 

i 1 


1 1 

1 1 


0.2 














1 V 


\ \ \ 


5 z 



Figure 9 

Precise solutions and iterations with decompositions with the cut-off. Here 
xiirn = 2 (the worst situation). Iterations are going near precise solution. 

With the growth of iteration number the lenght of coincidence between the 
iteration and the precise solution grows also rather rapidly. 
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